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Abstract. We consider classical and quantum one and two-dimensional systems with ladder 
OO ■ operators that satisfy generalized Heisenberg algebras. In the classical case, this construc- 

tion is related to the existence of closed trajectories. In particular, we apply these results to 
' (— I ', the infinite well and Morse potentials. We discuss how the degeneracies of the permutation 

^ i' symmetry of quantum two-dimensional systems can be explained using products of ladder 

operators. These products satisfy interesting commutation relations. The two-dimensional 
Morse quantum system is also related to a generalized two-dimensional Morse supersymmct- 
ric model. Arithmetical or accidental degeneracies of such system are shown to be associated 
to additional supersymmetry. 
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2 ; 1 Introduction 

Considering one-dimensional (ID) quantum and classical [2] Hamiltonians with polynomial 
ladder operators (i.e. polynomial in the momenta) satisfying a polynomial Heisenberg alge- 
■ bra [3113], we have recently been able to construct two-dimensional (2D) sup erintegr able systems 

^ , with separation of variables in cartesian coordinates, their integrals of motion and polynomial 

symmetry algebra. We have also discussed [5] how supersymmetric quantum mechanics [6] can 
be used to generate new 2D superintegrable systems. 

Sup erintegr able systems possess many properties that make them interesting from the point of 
view of physics and mathematics. Indeed, all bounded trajectories of maximally classical super- 
integrable systems are closed and the motion is periodic. Moreover, quantum 2D superintegrable 
systems have degenerate energy spectra explained by Lie algebras, infinite dimensional algebras 
or polynomial algebras. For a review of 2D superintegrable systems we refer the reader to [7]. 

Many ID quantum systems with nonlinear energy spectra such the infinite well, the Scarf and 
the Poschl-Teller systems have ladder operators that satisfy a generalized Heisenberg algebra 
(GHA) ilElIinillllliaiiaillliaiinillTlIISlIISlEO]. These systems also appear in the context 
of ID supersymmetric quantum mechanics (SUSYQM) [6]. It was also pointed out how some 
2D generalizations of the Morse potential are related to 2D SUSYQM |21j . 



*This paper is a contribution to the Proceedings of the Workshop "Supersymmetric Quantum Me- 
chanics and Spectral Design" (July 18-30, 2010, Benasque, Spain). The full collection is available at 
|http://www.emis.de/journals/SlGMA/SUSYQM2010.html| 
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At the classical level, it was recently discussed how ID systems like the Scarf, the infinite 
well, the Poschl-Teller and the Morse systems allow such ladder operators that satisfy a Poisson 
algebra that is the classical analog of a GHA [22 • An interesting paper [23] compared the ID 
classical and quantum Poschl-Teller potentials, their ladder operators and GHA. These ladder 
operators are directly related to time-dependent integrals of motion that give the motion in the 
phase space. In quantum mechanics ladder operators are well known and are used to generate, 
in particular, energy eigenstates, coherent and squeezed states [T8 l [2 ^ [25 l [26]. 

In this paper we are extending the constructions discussed in [H [2] to the infinite well and 
Morse classical and quantum systems. Such systems are integrable but not superintegrable. 
However some of their classical and quantum properties can be explained algebraically from 
their ladder operators as it is the case for superintegrable systems. 

In Section [21 we recall the trajectories of the ID classical infinite well and Morse potentials 
and we give some examples of closed trajectories for 2D such systems. In Section [S] we present 
general ladder operators and GHA for the ID classical systems and give them explicitly for our 
specific examples. We consider the extension to 2D of these systems and introduce some products 
of ladder operators which are not integrals of motion but satisfy interesting Poisson commutation 
relations. In fact, the commutation relations between these new operators can be interpreted 
as a condition to get closed trajectories in 2D. In Section [J] we recall the definitions of ladder 
operators and GHA for general ID quantum systems. We consider again the examples of the 
infinite well and the Morse systems [181 [211 [23 [221 [2Z] and we obtain the GHA. We also discuss, 
in Section [5] how using ladder operators, we can describe some degeneracies appearing from the 
permutation symmetry of isotropic 2D infinite well and Morse potentials. The interpretation 
of accidental or arithmetical degeneracies from an additional supersymmetry is discussed in the 
case of the 2D Morse potential. We end the paper with some conclusions. 



2 Classical trajectories 
2.1 Infinite well system 

The trajectories for ID and 2D infinite well systems are given in details in [28]. Here we just 
summarize the results to be able to use them later. In the ID case, we take the Hamiltonian: 



H, = ^ + V{x), V{x) 
2m 



0, \x\ < -, 

, , L 
oo, \x\ > —. 



(2.1) 



The trajectories are given by 



xit) 



to<t< ti, 



m 

Po \ ^ 

[t — t2u~l) + TT' *2i/-l <t< t2u, 

m 2 
m I 



(2n-l)mL 



Po 



(2.2) 



. The initial conditions are 



where n, G N* (i.e. are strictly positive integers) and t„ 

chosen as xq = a;(0) = and po = p(0) > at t = to = 0. The motion is periodic with period 
T = and frequency 



OJ 



vr 2E 




2mL L\m 

which is also written in terms of the total energy E of the system under consideration. 



(2.3) 
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For the 2D system, we take, in particular, the Hamiltonian: 



(2,4) 



with 



V{x) 



0, 



oo, 



|x|<-, 

, , L 
\x\>- 



V{y) 



oo. 



\y\ < 2 ' 
\y\ > —■ 



Trajectories are obtained directly from the ID case and periodic orbits occur when u^nx 



ylUy 



where Ux and ny E N* |28j . Using the equation (j2.3p adapted to 2D, we obtain the following 
relation 



Ex 



Ey. 



(2.5) 



We show in Fig. [T] a closed trajectory for the 2D infinite well taking the initial position as 
{xqiVo) = (0,0) and the initial momentum as po = {PQ,x-,Po,y) with po,a; > 0, po.j/ > 0. The 
parameters of the system are chosen as m = = 1. 



2.2 Morse system 

For the classical Morse system in ID the trajectories have been obtained in [29]. Let us recall 
some of the results. The Hamiltonian is given by 

^''%Vb(e-2/3^ -2e-^^). (2.6) 



Hx 



2m 



The trajectories can be easily obtained from the Newton equation. For negative energy E we 



get: 



1 — cos 6 cos{ujt + ^o) 



sm 



sm ( 



-E 



where is a constant depending on the initial conditions. In the case of positive energy, the 
classical motion is given by a similar expression involving hyperbolic functions. Let us mention 
that the initial position xq is different from zero while the initial momentum pQ is equal to zero. 
The motion for this ID system is periodic and can be viewed as the logarithm of a harmonic 
oscillator. Indeed, the frequency at energy E is given by 



uj = (3 



-2E 



m 



(2.7) 



We can easily extend this resolution to the 2D system given, in particular, by the Hamiltonian: 



H — Hx + Hy 



p2 p2 



2e-f'^)+Vo{e-^f'y -2e-^y). 



{21 



As for the case of the infinite well, closed trajectories and periodic orbits are obtained when 
^^xTT-x = ^yiT'y or, from equation (j2.7p . when the following condition is satisfied: 



n. 



■.■Sj -Ex = Uyy^-Ey, Ux^Uy^'M*. 



(2.9) 



We obtain here the analog of Lissajous figures for the anisotropic harmonic oscillator. One 
example is given in Fig. [2] where we have chosen the energies Ex and Ey higher than the 
minimum of the potential and such that equation (|2.9|) is satisfied. The parameters of the 
system are taken as m = /3 = 1 and Vq = 15. 
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Figure 1. A trajectory for Ej; — 2, Ey — i. Figure 2. A trajectory for E^ = —2, Ey 

3 Classical systems and generalized Heisenberg algebras 



We start by summarizing the factorization method in ID classical mechanics \22\ I23j . Let us 
assume that the system is given by the Hamiltonian H and admits ladder operators of the 
following form 

= ±if{x)P^ + g{x)VH + ^p{x) + ^(H), (3.1) 

where the real functions f , g, are supposed to depend on x only and is a real function of H 
(more precisely it is a power of ^fH). Let us also assume that the set {iJ, j4+,j4~} generates 
a generalized Heisenberg algebra given as 



{H,A^]p = ±i\{H)A^, 
{A'^,A-}p = -i^Ji{H), 
A^A~ = H --f{H), 



(3.2) 
(3.3) 
(3.4) 



where j{H), n{H), X{H) are power of VH. Note that the equation (j3.2p is a property that 
ladder operators must satisfy while the equation (j3.4p indicates a type of factorization of the 
Hamiltonian. In the case of bounded motions with negative energy, we will replace the square 
root y/H by ^/^. 

Let us mention that additional Poisson commutation relations are satisfied by the ladder 
operators of many classical systems \22\ [23] 

{{H,A^}p,A^}p = -2a{A^f (3.5) 

or equivalently 

{X{H),A^}p = ±2iaA^. (3.6) 
From equation (j3.2p . it is thus possible to construct two time-dependent integrals of motion 
Q± = A±e^*^(-^)*, (3.7) 
that indeed satisfy 



0. 
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The frequency of the bounded states has thus an algebraic origin and is given by X{E). The 
constant values of these integrals will be denoted by 

= c{E)e^''\ c{E) = ^E-j{E). (3.8) 

Indeed, the function c{E) is obtained from the equation (|3.4|) using the explicit form of the 
ladder operators given by the equation (j3.ip and the form of Q in equation (j3.7p 



= ±ifix)P^ + gix)VE + ij{x) + <I){E) = c(^)e±*(^o+^(^W. (3.9) 

The trajectories in the phase space {x{t),Px{t)) can thus be obtained algebraically from equa- 
tion (13. 9p with c{E) = (13. Sp . as we will show in the examples below. 

Now if we generalize the preceding considerations to a 2D system where the Hamiltonian 
allows the separation of variables in Cartesian coordinates 

H{x,y,Pa:,Py) = H^ + Hy, 

the system is integrable and possesses a second order integral of motion given by S = — Hy. 

Introducing the generalisation in 2D of the ladder operators (|3.ip and denoting them A^ 
and A^, we can form the following products 

= iAt)HA-y, /i^'') = {A-)HAty, 

where k and / G N*. They satisfy 

{H,lt'^}p = ±i{kXAH.) - lXy{Hy))I^^^'K (3.10) 
We have seen in [3] that for the special cases when \x{Hx) and Xy{Hy) are constant quantities 

(k I) 

and their ratio is rational, the functions ' are thus integrals of motion. In such cases, systems 
are superintegrable and all bounded trajectories are closed. More generally, the condition for 
which the right side of equation (j3.10p vanishes corresponds to the condition for the existence 
of closed trajectories (see equation (j2.5p or (j2.9p in our particular examples) 

kXxHEcc) - l\y{{Ey) = 0. 
The 2D extension of the additional relations ()3.6p implies the following constraints on /^''^ 
{{{//,/i'^'')},/i'='')},/J'')}^ = 0. (3.11) 

3.1 Classical infinite well 

The classical infinite well may be included in the preceding algebraic scheme {E > 0). Indeed, 
the ladder operators (with a slight modification) can be obtained from the quantum case |161 
[TSt [19] or as a limit of the classical Poschl-Teller system |221 [23] . The explicit form of the ladder 
operators for the infinite well depends on the boundary conditions. Indeed, using the same 
boundary conditions as in Section [21 we get 

A^ = T^=cos — P, + sin — /ff, (3.12) 
V2m \ L / \ L / 

which generate the generalized Heisenberg algebra given in equations ()3.2p . (13. 3p and ()3.4p with 
the following functions of H: 

\{H) = ii{H) = a^H , -f{H) = 0, a = — =. (3.13) 

Lv2m 
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2 

These ladder operators satisfy the equations (j3.6p (with a = The trajectories can thus be 
obtained algebraically from equations (|3.9p with c{E) = \fE. Indeed the real part leads to 

and thus, for the initial condition xq = and the phase choice = — f > we get 

x{t) = \ —t. 
V m 

We thus recover the equation (12. 2p . From the expression (12. 3p of the frequency, we indeed see 
that X{E) is equal to w. 

Let us now consider the following 2D classical infinite well as in equation (|2.4p . We obtain 
from equations (jS.lOp and ()3.13p 



{H,ir'}p = TMkVH^-ly^)ir- (3.15) 

These functions I^'^^ satisfy the equations (j3.1ip . We have thus shown that all closed trajectories 
have an algebraic origin and the equations (|3.15p vanish when the equation (12. Sp is satisfied. 

3.2 Classical Morse system 

The ID Morse potential given by equation (j2.6p possesses the following classical ladder opera- 
tors [22] (E < 0) 



We thus get the generalized Heisenberg algebra with the following functions of H: 

\{H) = eV^, ^[H) = e—j=, ^{H) = Vo + H + ^, e = -^. (3.16) 

Hy—H rl y Im 

2 

These ladder operators satisfy the equations (j3.6p (with a = — ^)- 

I 

We also have c{E) = y — -g — Vq. We thus get the trajectories algebraically and the fre- 
quency for the bound states are given by X{E) which is identical to equation ()2.7p . 

Let us now consider the 2D Morse potential given by equation (12. Sp . We obtain from equa- 
tions (IXTOD and (I3T6D 

{H,l'i^'^]p = ±ie{ky^,-ly^)I^^''\ (3.17) 

The right side of the equation (j3.17p vanishes again when equation (j2.9p is satisfied. All closed 
trajectories have an algebraic origin. 

4 Quantum systems with generalized Heisenberg algebras 

In ID quantum mechanics a generalized Heisenberg algebra [HI [9l [TO l [TT t [T2l [T3 l [T ^ [T5 l [T6t [T7 1 [TH l 
UHlllO] generated by ladder operators and Hamiltonian H is defined as the set {H, A'^ , A~} 
that satisfies the commutation relations 

[H, A+] = A+X{H) = 5{H)A+, (4.1) 
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[if, A-] = -X{H)A- = -A-5{H) (4.2) 

and 

[A+,A-]=j,{H)-^2{H), (4.3) 

with 

A+A- = ^i{H), A-A+ = -i2{H). (4.4) 

The functions X{H), S{H), ji{H), 72 (ii) will take special expressions depending on the 
physical system under consideration. 

It is well-known that, from equations (|4.ip or (|4.2p . we can obtain the energy spectrum 

algebraically. Indeed, let us consider the all set of eigenfunctions G 1} {I is a finite or 

infinite set of positive integers) such that 

H\,Pn) = En\4>n) (4.5) 

and 

A-\^n) = AARlV'n-i), A+\',Pn) = ^A:(n + 1) | . (4.6) 
From equations (j4.1j) . (14. 2p and (j4.6p . we get 

^{Efi) = '^(-E'n+l) = -E-n+l — En 

and, from equations (|4.4p and (j4.6p . 

7i(E„) = A;(?i), 72(^„) = /c(n + l). 

For the physical systems we are considering, the energy spectrum is at most quadratic. This 
implies that if we take 

En = ar? + bn + c, (4.7) 

we easily get 

A(£'„) = 5{En) + 2a = a{2n + 1) + 5. (4.8) 
Moreover, it is possible to show that, 

[[H,A^lA^] = 2a{A^f 
or, equivalently, 

[\(H),A^] = ±2aA^, (4.9) 



which is the quantum analog of the constraint (|3.5p or (|3.6p . 
Let us mention that if we take 

and assume that 

7i(i/)- 72(F) = -A(ii), (4.10) 
a 
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the ladder operators A^, together with the new generator Aq, generate a su{l, 1) algebra 

[A^,Ao]=±A^, [A+,A~] = -2Ao. 

When we act on the eigenfunctions IV'n), we also get 

[H,{A^)'^] = {E^±k-En)iA^)K (4.11) 

Now we consider again 2D systems as a sum of ID systems that have ladder operators 
generating the preceding algebraic structure. The Hamiltonian is written as: 

H{x,y,P,,Py) = H,{x,P,) + Hy{y,Py). (4.12) 

It admits a set of eigenfunctions labelled as {IV'n^.ny >,nx,ny £ 1} such that 

H\lpn^,ny) = En^^ny\lpn^,ny) with En^^Uy = + Eny 

As in the classical case, the separation of variables allows the existence of a second order 
integral of motion S = Hx — Hy. We can again take the products of the ladder operators as 
{k,leW) 

it'^ = {At)\A-)\ /i'^'') = {A-)\A-^J. (4.13) 

Let us impose the following constraints on the ladder operators (as in the ID case with 
equation (14. 9p ) 



[\{Hx),At] = ±2aAt [KHy),A^] = ±2aA^. 
From equations ()4.1ip and ()4.13p . we get 

[H, if^-')] = ±I^l''\k\x{Hx) - lXy{Hy) + k{k - l)a + /(/ + l)a). (4.14) 
Indeed, using equation (j4.1ip . we get, for example, 

[H, li'''^]\^n.,ny) = {{E^^+k - EnJ + (i?„,-; " i^n, ))/?'' Vn.,n, ) ■ 

Due to the fact that the total energy En takes the general form (14. 7p and \{En) is thus given 
by equation (j4.8p . we easily get 



{En^^k -EnJ + {Eny^i - Eny) = k{2anx + 6) - l{2any + h) + a{k'^ + l^) 

= k{Xx{EnJ - a) - l{Xy{Eny) -a) + a{k^ + f). 

This is the quantum equivalent of the equation (jS.lOp for systems with quadratic energy spec- 
trum. Similarly to the classical case, when \x{Hx) and \y{Hy) reduce to constants and their 

(k I) 

ratio is rational, the operators ' are thus integrals of motion and the Hamiltonian given by 
equation (14.121) is superintegrable [HE]. 

For our specific case of a quadratic spectrum, we have seen that the expression of A is given 
by equation (j4.8p and the commutators [H,!^'''^] are zero if a(2A;na; — 2/ny + /c^ + /^)+6(A; — /) = 0. 
In particular, for k = I and a 7^ 0, we get Uy = Ux + k. 

We will discuss, in the next section how, when \x{Hx) and \y{Hy) do not reduce to a constant, 
we can explain some of the degeneracies of the energy spectrum. 

Let us finally mention that we also obtain from equation (I4.14p the following commutation 
relations 

which are the quantum analogs of relations given by equations (jS.lip . 
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4.1 Infinite well potential 

We consider the quantum version of the ID infinite well given in equation (j2.ip . The Schrodinger 
equation and the corresponding energy spectrum are well-known and given by (using a as defined 
in equation ()3.13p ) 

H\ilJn) = Enlipn) , £'„ = , n = l,2,.... 

To simplify the following expression, we introduce the number operator N that satisfies A^IV'n) = 
n\^n) and is thus defined in terms of H as 

n = —Vh. 

an 

It has been shown that the ladder operators may be defined as [161 [TSl [19] 



2 f iL /"Kx\ /vrx 

A" = \ 1 + — { P^cos — - iVsin 



N \ TT \ L J \ L 



Let us mention the similarity with respect to the classical case and the equation p.l2p . These 
ladder operators satisfy the relations (14. 6p with 



k{n) = — 1. 

We can thus form a generalized Heisenberg algebra given by (j4.ip . (j4.2p and (|4.3p with 
\{H) = ahVH + ^ = ^{2N + 1), 

Note that the operator X{H) satisfies the relations of commutation given by equation (j4.9l 



with a = We also see that the relation (j4.10p is satisfied giving rise to the well-known 

existence of a su{l, 1) algebra associated to the infinite well system. 

4.2 Morse potential 

We consider the ID quantum version of the Hamiltonian (j2.6p . Introducing the parameters 



8mVo u-l 
we get the Schrodinger equation and the corresponding energy spectrum 



iJlV-n) =i^„,|^n>, ^^ = ^^[-^-'') = (4-15) 

where n = 0,1, . . . ,[p]. Indeed, the last admissible value of n is the integer part of p. Again the 
number operator N satisfies N\^pn) = nlipn) and is defined in terms of H as 

N = p- ^^/^. 
he 
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We introduce the following ladder operators [24:\ [25] 

= (v^)-' f (-^ft + f (2p - 2iV)(2p - 2JV - 1)) - (p + i) 



where 

^ (2p-Ar)(2p-2iV + 2) 
^ ^ (2p-iV + l)(2p-2iV)' 

They act on the eigenfunctions of the Morse potential as in equation (|4.6p where 
k{n) = n{2p — n). 



We can again form a generalized Heisenberg algebra given by equations (|3.2p . (j3.3p and (|37 
with 

A(i?) = /le^/^ - ^ = -^(2^ - 2iV - 1), 

71 (^) = - ^0), 12{H) = 71 (if) + j^X{H) 

and the operator X{H) satisfies the equation (I4.9P with a = ~ . Here again, we can construct 
a corresponding su(l, 1) algebra. 



5 Degeneracies of quantum systems 
with generalized Heisenberg algebras 

The study of the relation between the degeneracies and symmetries of multi-dimensional systems 
consisting in a sum of ID systems with nonlinear energy spectrum |30] appears to be a less 
studied subject than in the case of linear energy spectrum. This is a consequence of the fact 
that multi-dimensional systems that are sum of ID systems with linear spectrum possess the 
superintegr ability property [IllSlOlT] and appear in the classification of superintegrable systems. 
The most well-known of such systems is the anisotropic harmonic oscillator. 

Only the 2D anisotropic Q-oscillators was discussed in [20] and the permutation symmetry 
was obtained algebraically using operators {/+,/_, Is} (defined as series of ladder operators) 
generating a su{2) algebra and commuting with the Hamiltonian. This symmetry manifests 
itself only in doublet and singlet states. 

We will see how, from results of Section [4. 21 such type of operators may be considered as well 
for the cases of the isotropic 2D infinite well and Morse potentials. 

These last systems present also accidental or arithmetic degeneracies [30] that has not been 
explained algebraically. Some of the arithmetical degeneracies of the 2D infinite well were 
presented in [24] and references therein. We are proposing an interpretation for the case of the 
Morse potential based on the existence of supercharges of second order. In fact, this example 
is particularly interesting because this interpretation will also be valid for the super partner of 
the 2D Morse potential which is not separable in cartesian coordinates (not even in any either 
coordinates). 

For the 2D quantum systems we are considering, we have seen that they admit a quadratic 
energy spectrum of the form 

En^,ny = En, + = a{nl + n^) + b{nx + Uy) + 2c. 
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In particular, for the isotropic infinite well, we have a = -^-^ and b = c = while for the 

Morse system, we have a = c = =^ and b = The energies of these systems have similar 

algebraic structure and present the two types of degeneracies mentioned before. The first type 
appears when we make the change rix ^ Uy. It is associated to the permutation symmetry. The 
second type corresponds to the fact that two or more different sets {{ni^x,ni^y),i € N*} could 
give rise to the same energy. It is called accidental or arithmetic degeneracy |30] . 



5.1 Permutation degeneracies 

Let us consider the Fock space and the action of the ladder operators as given in equations (j4.5p 
and ()4.6p generalized for 2D. We use here the simplified notation \nx,ny) for the general energy 
eigenstates l^'n^^.n^)- We thus denote the operators J^^'^-* defined in ()4.13p by when k = I. 
We see that the operator takes a state with Ux = i and Uy = i + I to a state with = i + l 

and Uy = i and the operator I^^ takes a state with Ux = i + 1 and 7iy = i to a state with Ux = i 
and Uy = i + I. 

For the infinite well, we thus consider the following operators: 

1=0 1=1 ^ '' 1=0 l=\ 



/_ = /| = j:E^^N,z + o(M + ^i(/i'^) = EEiM + o(^+Mi, 

i=0 l=\ ^ ^ i=o 1=1 

and 

^ oo oo 
1=0 1=1 

They commute with the Hamiltonian of the system and describe the permutation symmetry 
in the energy spectrum of the 2D isotropic infinite well potential that manifests itself only in 
doublets and singlets. Indeed, we get 

I+\j,k) = \k,j), I.\j,k) = \k,j), h\j,k) = 0. 

For the Morse potential, the number of bound states is finite, we thus take only a finite sum 
for such kind of operators. We get 

[p] b]-« [p] \p]-i 

^+ = E E ) 1^' ^ + ^ + ^1 = E E 1^ + ^) ^ + ^1 ' 

1=0 1=1 ^ '' i=0 1=1 

\p] \p]-i [p] [p]-« 

= 4 = E E wTiv\^^ i + (M + ^1 (/i'^) = E E 1^' ^ + + Ml , 

i=0 1=1 ^ 1=0 1=1 



oo oo 



and 



^ [p] [p]-« 

-^3 = ^E E(i' + ^'^)^*+^'^i ~ \hi + i){hi + i\)- 



2 

i=0 1=1 



These products of ladder operators can thus be related to symmetries of the system as it is the 
case for superintegrable systems. 
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5.2 Arithmetical degeneracies 

We are considering the example of the 2D Morse system in order to give an interpretation of 
arithmetical degeneracies in terms of symmetry or more precisely supersymmetry of the system. 
Indeed, the existence of a super symmetric partner of the original system has been shown by 
loffe and collaborators (see [21] and references therein for more details and generalizations). We 
will thus start this subsection by introducing such a partner and the corresponding supercharges 
which realize the intertwining }21| . 

From equation (|4.15|) generalized for 2D, we get the energy spectrum of the Hamiltonian 
H = Hr, + Hy as 

En^,ny = + E^y = ^ ((p " n^f + (p " riyf), (5.1) 

with rixjUy = 0,1,2, ... ,[p]. The super partner of H has been obtained as 

/l2 /32 



H = H + 



m2sinh2ffal 



Indeed, introducing the supercharges Q^, which are differential operators of order 2 in the 
momenta, 

= -H, + Hy + D^, 
where are first order differential operators of the form 

we thus know that H and H are related by the intertwining relations: 

HQ+ = Q+H, HQ- = Q'H. (5.2) 

It means, in particular, that if ip is an eigenstate of H, Q'^ip is an eigenstate of H with the 
same eigenvalue. In fact, it has been proven that the new Hamiltonian H shares part of the 
energy spectrum of H. More precisely, the only normalizable eigenstates of H are given as (to 
simplify the developments, the set {nx,ny) has been renamed {n,m)) 

= Q+^i^ = {Em - En)^l„, + D+^>1^, (5.3) 

where 

"^It = ^«(rE)C(y) ±C(^)C(2/)), 

with ip^{x), V'md/) well-known eigenfunctions of the ID Morse system [18]. The func- 
tions (j5.3p are in fact antisymmetric in (n,m). The corresponding eigenvalues are given by 
En,m = (15. 1|) . for n,m = 0, 1, 2, . . . , [p] but now n ^ m and n 7^ m it 1 since ^^.n — ^-^^ 
^nn+i ~ ^- been shown [21] that we thus get the complete spectrum for H that is 

a subset of the spectrum of H. 

Let us summarize the properties involved in this context that will be useful for our interpre- 
tation of the arithmetical degeneracies. From ()5.2p . we get [Q~Q^,H] = and Q~Q'^'^n,m. = 

rn,m'^n,m with 

rn,m = (Em+i " En){Em-i " ^n) = {{m - nf - l) [{2p - m - nf - l) , (5.4) 
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where r„,n+i = '''n+i,n- In fact, Q Q'^ is a fourth order differential operator acting on the 
eigenstates ^n,m of H that could be written as 

Q-Q+ = {H^ - Hyf + 2H + I. (5.5) 

Moreover, we get [Q'^Q'^^H] = and Q'^Q^^^j^ = rn^m^nm- The operator Q'^Q~ is again 
a fourth order differential operator acting on the eigenstates of ^■ 

The interesting consequence is that we have found a superalgebra constructed from both 
Hamiltonians and the corresponding supercharges as follows. We introduce the generators 

that satisfy 

[n,Q^] = o, {Q+,Q-} = 7^, 

where TZ commutes with Ti. Thus the set {H, TZ, Q^, Q~} closes a superalgebra. Note that TZ is 
a fourth order differential operator which will be useful to explain the arithmetical degeneracies 
of the spectra of both H and H. 

We are now ready to proceed to the analysis of arithmetical degeneracies of the 2D Morse 
systems. Let us consider, for example, the Hamiltonian H and assume that it admits a double 
arithmetical degeneracy in the energy spectrum for the couples (ni,?Tii) and (ra2,m2). It means 
that -Em, mi = -E'n2,m2 (avoiding the permutation symmetry and the identical relation). From 
the supersymmetric context, we have introduced a new constant of motion TZ with eigenvalues 

' n,m — 

n . We are thus showing that rni,mi 7^ i^n2,m2- 
Let us first set p — rii = ki and p — nii = U. Since En^^m^ = En2,m2i we get 

If we express ki in terms of the other quantities, we see that the difference ^i — '''n2,m2 
becomes 

'"ni,mi ~ fn2,m2 — ~ ^2) (^1 ~ ^2) 

= 4(m2 - mi){n2 - mi){2p - (m2 + mi))(2p - (n2 + "ii)). 
To get '"ni,mi ~ fn2,m2 — 0; wc have 4 possibilities: 

1) mi = m2 and this leads to Em^mi = En2,mi <^=^ ni = n2 (identity, excluded); 

2) 772-1 = n'2 and this leads to -Em, mi = -E'mi,m2 '^i = ^2 (symmetry permutation, 
excluded) ; 

3) 777-1+7722 = 2p and this leads to r?7i = r?72 = p since 7?7,j < p and En-^^^p = En2,p <^=^ ni = 722 
(identity, excluded); 

4) r?7i + 772 = 2p and this leads to r?7i = 722 = ^' and En-^^p = -£^p,m2 *^=^ ^1 = "72 (symmetry 
permutation excluded). 

So we conclude that r„,^,mi 7^ fn2,m2 and the arithmetical degeneracies are explained by the 
existence of TZ since the preceding proof may be extended to multiple degeneracies. 
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6 Conclusion 

In earlier works [H [21 [3] , we have presented some ways of getting integrals of motion and poly- 
nomial symmetry algebras in classical and quantum mechanics for multi-dimensional superinte- 
grable systems. These quantities are obtained from ID systems with polynomial ladder operators 
satisfying polynomial Heisenberg algebras. The integrals are products of ladder operators. 

It was also observed in context of superintegrable systems that higher order ladder operators 
and thus higher order integrals of motion can be written in terms of first order and second order 
supercharges [U [U [31]. In particular, we constructed higher integrals of motion for an infinite 
family of superintegrable systems and involving the fifth Painleve transcendent using products 
of second order supercharges [31]. However, all these Hamiltonians were separable in cartesian 
coordinates. Let us mention that it has been shown very recently that a new non separable 
superintegrable system admits one third and one fourth order integrals of motion |32j . 

In this paper, we have considered systems with non polynomial ladder operators. The 
method does not produce superintegrable systems and this is a limitation to the construc- 
tion of [H [21 [3]. However, 2D integrable systems with non polynomial ladder operators can 
be constructed. In the classical case, we have obtained algebraically the trajectories and the 
condition for the existence of closed trajectories. In the quantum case, we have explained the 
degeneracies of the energy spectrum which involve, in particular, supersymmetric methods. 
We have applied these results for the classical and quantum infinite well and Morse poten- 
tials. 

Finally, let us observe that the new integral of motion TZ, which is a fourth order differential 
operator, is given as 



n 



Q+Q~ 




and acts on the superspace of eigenfunctions {(0, (^^^^i ^nmY' ^ n > m; n,m = 0,1, 2, 

...,\p]}. We see that Q~Q^ = (jS.Sh is quadratic in and Hy. The operator Q'^Q' has 
a similar expression but with additional terms. Indeed, we can show that 

Q+Q- = {H, - Hyf + 2H + I+ [if, - Hy, D+ - D']. 

We could say that we have a nonlinear supersymmetry but, in fact, these last operators are not 
quadratic in H ox H. 

We see that the supercharges are related to = Hx — Hy , the integral of motion obtained 
somewhat trivially for Hamiltonians admitting the separation of variables in cartesian coordi- 
nates. The fourth order integral of motion IZ associated to the Morse Hamiltonian in 2D and 
related to the "square" of Q is very interesting since one of the Hamiltonian {H) does not allow 
the separation of variables. 
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